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ABSTRACT 


In August of 2006 the Adaptive Sampling and Predietion (ASAP) experiment was 
eondueted near the northern Monterey Bay. Multiple assets including aircraft, 
autonomous vehicles, moorings, and numerical models were used to gain a better 
understanding of three-dimensional upwelling centers. Data were collected at two 
separate mooring locations using Acoustic Doppler Current Profilers (ADCPs) during the 
experiment. The focus of this thesis is to determine the effects of local wind forcing on 
the ocean circulation and provide a comparison between the data collected at the mooring 
locations and numerical predictions for the region. Upwelling and relaxation events are 
used as the basis for understanding the local wind forcing. Upwelling typically results in 
equatorward flow while relaxation events typically result in poleward flow. Several 
different types of analyses were used to determine the effects of the local wind forcing. A 
visual analysis was performed with stick vector plots and component plots of the rotated 
time series that compared the wind with the data from the water column. Two methods of 
cross correlation, component correlations and vector correlations, were exploited as well 
as a spectral analysis of the wind and ADCP data. Finally the coherence and phase 
between the wind and currents were examined. Based on the analysis it became evident 
that the currents were forced by both wind and non-local events such as eddies, 
meanders, and the large-scale alongshelf pressure gradient. 

Associated with the ASAP experiment, the Harvard Ocean Prediction System 
(HOPS), the Regional Ocean Modeling System (ROMS), and the Navy Coastal Ocean 
Model (NCOM) provided nowcasts that were compared with the mooring data to 
determine their accuracy and precision. Overall, in the beginning of August the models 
provided reasonable representations of the flow patterns at the mooring locations. The 
prediction error increased towards the end of August which was possibly related to data 
assimilation techniques and more non-local forcing at that time. The military application 
of this thesis is that accurate current prediction by ocean models will benefit amphibious 
operations, special warfare operations, and mine warfare in the littoral zone. 
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I. 


INTRODUCTION 


A. BACKGROUND 

Due to its proximity to both the continental shelf and the Monterey Bay 
submarine canyon, the Monterey Bay provides a unique area for oceanographic research. 
The waters off the California coast are also part of the California Current System, a 
prominent eastern boundary current. Eastern boundary currents are generally linked to 
major upwelling systems with circulation patterns similar to those observed in the 
vicinity of the Monterey Bay. Wind forcing experiments have previously been conducted 
at areas north of the San Francisco Bay because the continental shelf there has 
comparatively uncomplicated bottom topography, and changes due to the large-scale 
wind stress can be identified (Winant et ah, 1987). 

In August of 2006 the Adaptive Sampling and Prediction (ASAP) experiment was 
launched in northern Monterey Bay. The major goals of the ASAP experiment were to 
use multiple observational assets and data assimilation in an attempt to gain the best 
insight into model predictive skill. First, autonomous vehicles including two types of 
gilders were deployed in patterns surrounding the Point Ano Nuevo upwelling center. 
The mobility of the gliders makes them efficient sampling resources since they can 
observe the constantly changing ocean environment. These unmanned sampling devices 
can provide near-real time data as they sample the ocean. As the data are analyzed 
onshore, the paths of the gliders can be changed to sample interesting oceanic elements 
such as eddies and thermal fronts, a process known as adaptive sampling (Princeton 
ASAP Program, 2006). The real-time data were then assimilated into three separate 
numerical models focused on the Monterey Bay region. The overall goal is to have a 
system of data/model assimilations that produce predictions of oceanographic conditions 
with the smallest error. 

This thesis project is designed to analyze data recorded by two Acoustic Doppler 
Current Profilers (ADCP) moored in an across-shelf configuration as part of the ASAP 
experiment. The ADCPs served as points of comparison between the moorings, 
autonomous vehicles, and model outputs. The data from these profders were fdtered. 
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rotated, plotted, and eompared. The statistical analysis of the data helps to determine the 
effects of upwelling and relaxation events on the circulation of the water column north of 
the Monterey Bay. In particular, there is an episodic upwelling of cold water that 
historically appears there during the spring and early summer time frame. The data will 
also be compared with model nowcasts and forecasts from three ocean models associated 
with the ASAP experiment. 

B. FUNDAMENTALS OF SHELF CIRCULATION 

The area north of the Monterey Bay is a convenient area to conduct field work 
pertaining to wind forcing, since upwelling and relaxation events appear to be responsible 
for a good percentage of the coastal circulation in this region. The surface circulation in 
this area tends to be equatorward in the summer with evidence of eddies and meanders, 
tidal variabilities, and the local sea breeze (Ramp et ah, 1997; 2005). Through data 
collected by the ADCPs and the model/data comparisons, research is conducted to 
observe the effects of wind forcing on the speed and direction of ocean currents. This 
project attempts to describe key aspects of ocean circulation mainly focusing on coastal 
upwelling and relaxation cycles. The upwelling process is important because it stimulates 
mixing, brings nutrients and dissolved chemicals to the surface layers, and contributes to 
surface water temperatures and productivity. 

I. Wind Forced Circulation 

The earliest contribution to understanding the wind forced circulation was made 
by oceanographer V. W. Ekman (Ekman, 1905). Within the surface boundary layer 
surface currents are directed 45° to the right of the wind. This forcing is propagated 
downward throughout the water column until the current has been rotated 180° clockwise 
from its starting direction. The decreasing magnitude of the downward rotating current is 
known as the Ekman Spiral. The depth at which the current direction is 180° from the 
surface vector is referred to as the Ekman depth. The Ekman transport is the integrated 
transport over the Ekman depth, and is directed 90° to the right of the wind. As a result 
the equatorward winds off the California coast transport surface water offshore during 
upwelling events (Eigure 1). The offshore transport is balanced by onshore flow in the 
near-bottom layers. 
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Effects of wind forcing have been observed within the Monterey Bay and its 
vicinity. Strong upwelling favorable equatorward wind events lasting periods of two to 
three weeks are common during the upwelling season. These periods of strong winds 
typically lead to upwelling conditions and create a cyclonic circulation throughout the 
bay. Southward flow is found across the mouth of the bay, northwestward flow is found 
across the northern shelf, and the circulation is completed with eastward flow on the 
southern shelf (Ramp et ah, 2005). The circulation becomes more chaotic and less 
defined during relaxations in the winds. Additionally, the California undercurrent flows 
poleward across the mouth of the Monterey Bay during these relaxations (Ramp et ah, 
2005). 

The wind stress is typically effective in directly driving currents on the 
continental shelf due to the divergence of the Ekman transport at the coast and the 
shallow topography. The amount of correlation between the wind stress and circulation 
depends on a number of factors including the presence of a coastline, the bottom 
topography, and the wind near the coast as well as the presence of other significant 
oceanic features such as eddies or pressure gradients. 

2. Upwelling Events 

In general an upwelling event is characterized by colder bottom waters replacing 
warmer surface waters closer to shore, a process that occurs as the surface layer is forced 
away from the coast by the surface wind stress (Graham and Eargier, 1997; Hickey, 
1998) (Eigure 2). In a relaxation event there is poleward or no wind. Wind relaxations 
can usually be observed in the alongshelf component of the current flow (Send et al, 
1987). 

Along a straight coast, this two dimensional process can be described in terms of 
a vertical and across-shelf direction. Offshore Ekman transport causes the isopyncals to 
slope upwards toward the coast and generates an equatorward alongshelf jet. The 
upwards slope is associated with a mean vertical shear in the alongshelf current flow 
(Huyer, 1984). Within the California Current System, the upwelled waters are saltier as 
well as colder. Mixing and vertical advection can also occur creating a homogenous/less 
stratified environment near the coast. This newly defined thermal structure also affects 
the equatorward direction of the current. 
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During a coastal upwelling situation the sea level at the coast will decrease. This 
change in sea level, reflected in the southward flow of the alongshelf surface currents, 
ereates a barotropic pressure gradient directed eastward towards the coast (Huyer, 1984). 
At this point the eurrent flow near the coast is consistently downwind and not following 
the direetion of the Ekman spiral henee the transport is less than 45° to the right of the 
wind (Ekman, 1905). The total eurrent is therefore the sum of the equatorward 
geostrophic flow, the baroelinic jet, and the Ekman flow. As depth increases offshore, the 
alongshelf eurrents often reverse to form a poleward undercurrent flowing in a direetion 
opposite of the wind foreing (Kosro, 1987; Huyer, 1984). This circulation model is a 
typical result of coastal ocean dynamics. 

Many field experiments focusing on shelf circulation along a straight coast line 
have been eompleted in the vicinity of the northern California eontinental shelf and slope. 
In 1981 and 1982 the Coastal Oeean Dynamies Experiment (CODE) suecessfully 
reeorded two four-month data sets through two upwelling seasons with an emphasis on 
summer upwelling conditions. In relation to wind forcing, the CODE study found that the 
alongshelf wind is eoherent with the eurrents and appears to have a direet effect on the 
direction of current flow and temperature change (Send et al, 1987). Conducted between 
Pt. Arena and Pt. Reyes, the CODE study took place in an area with the strongest 
upwelling favorable winds in the California Current System (Huyer, 1984). The CODE 
study discerned that the combination of the strengthening of the North Pacific high and 
the resulting colder surfaee temperatures of approximately 9 °C announced the beginning 
of the upwelling season (Winant et ah, 1987). 

A variable eoastline geometry with headlands results in a three-dimensional 
upwelling cireulation. Vertical, across- and alongshelf directions are all considered in a 
three-dimensional model. In a three-dimensional cyele, the general meehanics of 
upwelling remain the same as compared to a two-dimensional cyele. The difference 
between the two is determined by the location of the colder upwelled water as a funetion 
of coastline topography. Reeall in a two-dimensional upwelling eycle the eolder water 
will be more or less evenly distributed along the straight eoastline. In a three-dimensional 
situation with headlands, the upwelled water is concentrated in a general location off of 

the headland, often referred to as an upwelling center. 
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Upwelling centers are regions of loeally intense upwelling, which have been 
observed at several loeations along the west eoast of the United States. These regions 
have a eharaeteristie “tongue” of eold water extending from the eoast to several hundred 
kilometers offshore (Kosro, 1987). A deerease in surfaee temperature results in 
diminished vertieal temperature, and eorresponding density gradients over the continental 
shelf (Winant et ah, 1987). Observing and understanding the properties of the water 
column will determine how the circulation in an upwelling situation is altered by outside 
forcing. 

There are two major upwelling eenters in the Monterey Bay region, loeated to the 
south at Point Sur (Breakers and Mooers, 1986) and to the north at Point Ano Nuevo. 
While Point Sur is a larger, more energetic upwelling center, this experiment focused on 
the upwelling eenter at Point Ano Nuevo due to its closer location to Santa Cruz and 
Moss Landing. The cold waters from the Point Ano Nuevo eenter have been observed to 
flow south and circulate within the Monterey Bay region, thus impacting local upwelling 
and temperature gradients. Speeifieally in the Monterey Bay region eold water first 
appears in the north but eventually spreads southward as a result of the wind (Rosenfeld 
et ah, 1994; Ramp et ah, 2005). 

The Monterey Bay region experienees eoastal upwelling for approximately six 
months from Mareh through August. The northern California seasonal upwelling begins 
following the annual spring transition that oeeurs when the large-scale wind beeomes 
equatorward during the spring (Hiekey, 1998). The six-month upwelling season ean be 
elassified into three different phases. In March when the upwelling season begins the 
Monterey Bay is eharaeterized by surfaee waters near 13 °C, a slight indication of surface 
warming, and low salinities (Graham and Largier, 1997). Colder waters (<10.5 °C) and 
higher salinity values (>33.7 psu) are indieations of the peak of the upwelling cyele 
occurring during April through June (Rosenfeld et ah, 1994). Towards the end of the 
season, July through September, warmer temperatures (13-16 “C) are observed in the Bay 
(Graham and Largier, 1997). 

3. Wind Relaxation Events 

In the northern California shelf region, relaxation events occur throughout the 

spring-summer upwelling season. It is intuitive to conclude that relaxation events are the 
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reverse of upwelling events. If this were an accurate situation then when the upwelling 
winds stop the warmer water that was dislocated offshore would move back to its initial 
position closer to the coast and the sea level would increase. This, however, does not 
always happen. 

During a wind relaxation event, warmer less saline water appears along the coast 
(Ramp et ah, 2005; Rosenfeld et ah, 1994). Even though the dominant process is not fully 
understood the appearance of warmer coastal waters after a relaxation event can be 
related to alongshelf advection or surface heat fluxes. As a result of this change, the 
isotherms and isohalines that were sloped upwards to the coast in the upwelling event 
adjust to a more level state. In this case across-shelf advection and flow do not contribute 
to the increased warm waters observed during relaxation events (Send et ah, 1987). 

In the Monterey Bay region, relaxation events are observed more frequently and 
have longer time periods towards the later part of the season (July and August). 
Prolonged relaxation events tend to result in shoreward movement of the California 
Current into the Monterey Bay causing low salinity values and warmer water 
temperatures varying around 15 °C (Graham and Largier, 1997). Previous studies of the 
Monterey Bay region have also observed rapid onshore advection of warmer waters 
combined with surface warming as far inshore as the beginning of the Monterey 
Submarine Canyon (Rosenfeld et ah, 1994). 

C. OVERVIEW OF MODEL SPECIFICATIONS 

Another important aspect of the ASAP experiment was to evaluate the importance 
of data assimilation to numerical predictions. The three models used during the ASAP 
study were the Harvard Ocean Prediction System (HOPS) (Robinson, 1999), the 
Regional Ocean Modeling System (ROMS) (Schepetkin and McWilliams, 2004), and the 
Navy Coastal Ocean Model (NCOM) (Shulman et ah, 2002). A comparative table of the 
numerical models is included for reference (Table 1). These nested mesoscale models 
were adapted to the area of interest off Point Ano Nuevo just north of the Monterey Bay. 
Even though the three numerical models are designed to predict the same variables, there 
are many factors that can cause differences in the outputs. Each model has a different 

domain and set of initial conditions used to set up the model. These domains are designed 
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to encapsulate the oceanic/atmospheric conditions of the particular area of interest. One¬ 
way nested models, such as the NCOM and ROMS models, have larger domains that 
provide boundary conditions for the smaller areas of interest (Figure 3). The ROMS 
model uses a triply-nested domain which covers the west coast, central California coast, 
and the Monterey Bay at 15 km, 5 km and 1.6 km resolution respectively (Jet Propulsion 
Lab, 2006) (Figure 3). Whereas one-way nested models only communicate from the 
larger domain to the smaller area of interest, two-way nested models can share boundary 
influences both ways between nests. The HOPS model in particular is two-way nested in 
a local scale model whose largest domain only covers an area between the Monterey Bay 
and San Francisco Bay. There were two separate HOPS runs used for comparison. The 
HOPS Solo model has 1.5 km resolution and the HOPS Nested model is a two-way nest 
with 1.5-0.5 km resolution. The two HOPS models also use different bottom friction 
parameters. 

Most oceanographic models use one of three coordinate systems to resolve the 
water column. Sigma coordinate models use vertical coordinates that follow bottom 
bathymetry. The sigma coordinate method is more effective in complex topographic 
situations, coastal regions and over the continental shelf. To most effectively map the 
water column, sigma levels are typically closely spaced near the surface and/or bottom 
(Figure 4). This separation scheme allows for improved resolution of the boundary layers 
(Chassignet et ah, 2002). The HOPS model uses 30 sigma levels to resolve the water 
column from surface to bottom. They are most effective in the surface mixed layer 
because they are capable of providing the high resolution needed to represent 3D 
turbulence (Chassignet et ah, 2002). The z coordinates use a predetermined number of 
depth levels to resolve the water column and do not provide accurate resolution in areas 
with sloping topography (Figure 5). There is also another method, which involves a 
hybrid system between the sigma and z levels to obtain the most effective representation 
of the water column. These models use sigma coordinates in the upper ocean and z 
coordinates in the deep ocean and away from abrupt topography (Figure 6). In this study 
all three models used the sigma coordinate system. The ROMS model uses 32 sigma 
levels and the NCOM model uses 30 sigma levels. 
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Mixing is one of the most eritieal modeling aspeets to consider because it 
provides significant differences between model outputs. It is necessary for a model to use 
an accurate range of mixing parameters. Too little mixing could cause a model to crash, 
while a large mixing parameter could produce model outputs that are too smooth. Due to 
the nature of this study and the focus on upwelling and relaxation events, mixing plays a 
large role in the model output. The NCOM model used the most common mixing 
scheme- the Mellor-Yamada 2.5 scheme, while the ROMS model used the K-Profile 
Parameterization. The HOPS model used a two-mode scheme involving a top layer that 
incorporates wind stress and Ekman depth and a bottom layer that is based on the local 
gradient Richardson number. 

In addition to mixing, models must also have a starting point. Each model has its 
own set of initial conditions for the predicted/output variables. The HOPS model used 
recent conditions from the local area to form initial conditions, while the ROMS and 
NCOM models used climatology that was forced by the Coupled Ocean Atmosphere 
Mesoscale Prediction System (COAMPS™). Generally, the background climatic 
information is entered into the models during the spin-up process. Spin-up describes the 
amount of time needed for a model to reach a state of statistical balance under forcing. 
This process builds a historical climatic reference of the region of interest and the 
surrounding area. Depending on the purpose and goals of the model the spin-up process 
can vary in time, the amount of included historical information, and the type of 
climatologic database. The ROMS model had a significantly longer spin up time, starting 
in July of 2003, than the NCOM and HOPS models started in July 2006. All three 
models used COAMPS data to define air-sea heat flux parameters as well as surface wind 
stress parameters. The ROMS model used the COAMPS data for wind information and 
used the Earge and Pond drag coefficient to compute the wind stress (Earge and Pond, 
1981). 

In an effort to provide the most accurate prediction, there are different data 

assimilation techniques utilized by models. The fundamental idea of data assimilation is 

to incorporate real data into the simulated ocean system that is predicted by the model. 

This is achieved by identifying the correlation scales, determining the correct amount of 

influence the real data retains, and finally applying the accurate amount of change to the 
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model output. Data assimilation is a delicate process since adjusting one predicted 
parameter in accordance with the real data will have a ripple effect on the remaining 
variables in the model output. 

The three models in this study apply three separate methods of data assimilation 
within continuous assimilation schemes. An assimilation window describes the regularity 
of the data assimilation. The HOPS and NCOM models used twelve-hour updates while 
the ROMS model used six-hour updates. All three models assimilated temperature and 
salinity from the gliders and sea surface temperatures from the aircraft. ROMS 
additionally assimilated satellite sea surface temperatures and temperate and salinity from 
MBARI moorings Ml and M2. The NCOM model did not apply satellite data but used 
temperature and salinity profiles along with sea surface temperatures from the Modular 
Ocean Data Assimilation System (MOD AS) (Summary of Three Ocean Models, 2006). 
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II. DATA COLLECTION 


A. ACOUSTIC DOPPLER CURRENT PROFILERS 

Current information for the 2006 eireulation experiment was reeorded using two 
bottom-mounted Teledyne RD Instruments Sentinel Workhorse 300 kHz ADCPs. 
Loeated at 36° 55.336'N, 122° 07.344'W, ADCP 1 was plaeed in shallower water of 54 
meters. ADCP 2 was loeated at 36° 53.683'N, 122° 11.244'W at a depth of 92 meters 
(Figure 7). The moorings were loeated 6.55 km apart. The wind time series was reeorded 
by buoy M2 loeated in the Monterey Bay and maintained by the Monterey Bay Aquarium 
Researeh Institute (MBARI). Buoy M2 rather than MI was used for the ASAP 
experiment as its loeation farther offshore was more representative of the winds at ADCP 
I and 2. Bottom-mounted moorings are advantageous beeause they have no extraneous 
movement, and allow faster sampling and less averaging than a surfaee mooring. The 
ADCP is eapable of reeording the data in four-meter bins up to the surfaee of the water; 
however, the sidelobes refleet from the surfaee and eontaminate the first three bins. The 
ADCP has a eorrelation magnitude feature that allows it to determine when the reeorded 
data is good or not. 

The ADCP uses a multiple beam system to obtain an aceurate observation of 
eurrent veloeity. The eurrent profilers are designed to sample the water eolumn every 
seeond from four separate beams on a 20° slant. Sinee only three beams are required to 
evaluate veloeity, the fourth beam serves as a way to assess the quality of the data 
(Teledyne RD Instruments, 2007). Three beams resolve the x,y,z direetions while the 
fourth beam allows eomputation of the error veloeity, the differenee between two vertieal 
veloeity estimates. This fourth beam is designed to assess aeeuraey, reliability, and 
quality of the data (Teledyne RD Instruments, 2007). 

The ADCPs were used for data eolleetion over a one-month period. Data were 
eontinuously eolleeted and an ensemble was ereated as the one-seeond samples were 
averaged over a one-minute period. The ADCPs were programmed to send baek the most 
reeent estimate via satellite link, every 10 minutes. 
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B. COORDINATE ROTATION 

To ease analysis, the prineipal axis teehnique was used to rotate the raw wind and 
ADCP data. Rotating the data separates the u and v eomponents into aeross- and 
alongshelf eomponents. This method of rotation ensures that a minimum of the 
alongshelf flow will be translated into the aeross-shelf flow (Kosro, 1987). The prineipal 
axis rotation agreed with the loeal orientation of the bathymetry and yielded angles of 58° 
for ADCP 1 and 48° for ADCP 2. 

C. DATA FILTERING 

The raw reeorded data from both ADCP mooring sites were truneated and 
eheeked for errors. Onee the data were transeribed into a form useable by MATLAB, 
they were filtered on three different time seales to inelude energy with periods less than 
five hours for internal waves, periods between five and thirty-three hours for tidal and 
inertial motions, and greater than thirty-three hours for wind foreing, mesoeale events, 
and mean currents. These different filtering bands allow different ocean phenomena to be 
observed in the plotted outputs. Time series plots were used extensively throughout this 
research for data analysis. Additionally, the filtered and rotated aeross-shelf and 
alongshelf component plots and color contours were also generated using MATLAB. 
These plots aid in interpreting statistical computations of the data. Autospectra were used 
to show the distribution of energy as a function of frequency and coherence and phase 
were computed between wind and currents as well as between the time series from the 
two ADCP moorings. 
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III. DATA ANALYSIS 


A. MEAN FLOW 

In order to achieve an accurate understanding of the circulation, it is necessary to 
analyze different depths of the water column. Table 2 outlines the depths used for the two 
moorings. As previously mentioned, the near-surface depth is not directly at the surface 
due to sidelobe interference. The near-bottom depth for ADCP 2 represents the deepest 
depth provided by all four models. 

The mean currents at ADCP 1 were all poleward. The surface flow (2.05 cm/s) 
was less than the subsurface flow (4.84, 4.19 cm/s) due to the opposing equatorward 
wind stress. The across-shelf means were small, less than 1 cm\s. At ADCP 2 the surface 
mean was equatorward (-6.68 cm/s), the bottom was poleward (1.69 cm/s) and the mid¬ 
depth mean was near zero. The across-shelf means were similar to ADCP 1, less than 1 
cm/s. The standard deviations at ADCP 2 were about double ADCP 1, indicating much 
larger current fluctuations there (Table 2). 

B. VISUAL ANALYSIS 

The velocity contours show significant flow variations within the water column. 
Three poleward flow events can be clearly observed. ADCP 2 showed strong poleward 
flow starting on 03 August and lasting to 08 August (Figure 8). This four-day event 
resulted in strongest flow between 05 and 06 August and occurred partially during a wind 
relaxation event. Poleward current flow may also be associated with alongshelf pressure 
gradients and mesoscale variabilities. The poleward event starts earlier than the wind 
relaxation event hence this strong forcing could be a coincidence or a combination of the 
wind relaxation and other oceanic processes. Additionally, there were two visibly weaker 
poleward flow events that took place later in the month during equatorward winds. These 
events appear to be more affected by remote forcing versus local forcing 
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The first poleward flow at ADCP 2 appears to lead a weaker flow event at ADCP 
1 by about one day (Figure 9). This is indieative of an oeeanic process that could be 
moving onshore such as an eddy or front. Additionally, shorter spatial and temporal 
scales were observed in across-shelf components at both mooring locations indicating 
evidence of wind forced currents in the surface layer. 

Before analyzing the directional effects of the winds and currents it is necessary 
to discuss the differences in current magnitude at the individual mooring locations. The 
upwelling winds are on the order of 5 m/s and occur at the same time as currents on the 
order of 5-10 cm/s at ADCP 1 and in excess of 15 cm/s at ADCP 2. There are several 
possible explanations to consider concerning the stronger current flow at ADCP 2. Firstly 
due to location, ADCP 2 was located closer to the center of the upwelling jet where 
ADCP 1 was located farther in s hore of the jet. A second possible reason for increased 
flow rates is related to the offshore increase in the wind stress, causing more wind stress 
at ADCP 2 than at ADCP 1. This spatial variation was observed in the aircraft data, for 
instance on 01 August 2006 (Figure 10). 

The stick vector plots clearly show the strong poleward flow from 03 August to 
08 August (Figures 11 and 12). After comparing with the wind events, it becomes clear 
that a wind relaxation occurred during a shorter time frame than the poleward flow. The 
wind relaxation is also consistent with the aircraft data taken during the experiment. 
Small or poleward wind velocities are evident on 03-04 August (Figure 13). It is evident 
at both mooring locations that the poleward current leads the relaxation event. As the 
time series progresses, both mooring sites show a transition day on 07 August where 
upwelling favorable winds start in the north and the flow of the poleward currents start to 
decelerate. Again this transition was also apparent on the aircraft data. Figure 13c shows 
the beginning of the equatorward wind shift. By 08 August, the aircraft data indicates 
strong equatorward upwelling favorable winds consistent with the stickplots (Figure 
13d). Around 09 August the currents become equatorward thus following the southward 
wind events indicative of an upwelling situation. 

A wind reversal was observed around 14 August but there was no significant 
visible change in the currents at any depths at either mooring. Lack of reversals in the 
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water column can be related to the required response time between currents and wind. It 
is observed that the current weakens only to increase again around 16 August in relation 
to the upwelling wind event. Equatorward flow was observed at both mooring sites and 
appeared to be stronger in the surface currents than the mid-depth or near-bottom currents 
for both ADCP 1 and 2 (Figures 11 and 12). A weak relaxation event on 18 August was 
followed by a poleward shift of the current. Despite strong upwelling winds, this 
poleward shift appears to dominate the next ten days at both moorings. It is interesting to 
note that the upwelling winds are the same strength (5 m/s) as they were on 08 August 
when they were accompanied by an equatorward flow. As observed at both mooring 
locations, the poleward flow from 18 August to 28 August was not as strong as previous 
flows. Late in the summer, there is often a poleward flow that is associated with an 
alongshelf pressure gradient force. 

The wind time series showed an upwelling event that began on 26 August and 
continued to the end of the time series. At ADCP 1 the surface currents showed weak 
equatorward flow that was less than 5 cm/s (Figure 11). The mid-depth and near-bottom 
currents showed no immediate change from poleward to equatorward flow until 30 
August. After the currents switched, the equatorward flow was weak, roughly 2 cm/s. 
There is also some previous evidence that the poleward pressure gradient force increases 
with the season (Ramp and Abbott, 1998). At ADCP 2, the flow at all three depths 
clearly became equatorward and followed the upwelling favorable wind (Figure 12). 
Overall, it appears that there was a stronger response at ADCP 2 to the wind forcing than 
ADCP 1. 

C. CROSS CORRELATIONS BETWEEN MOORING LOCATIONS 

Vector and scalar correlation values were calculated between ADCP 1 and ADCP 
2. Individual comparisons of the time series plots show that there are higher correlations 
associated with the alongshelf component than with the across-shelf component. In the 
alongshelf component, correlation increased proportionally with depth. Table 3 displays 
the cross correlation coefficients for the maximum across-shelf and alongshelf values as 
well as the respective significance levels for the three indicated depths. There is strong 
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correlation in the alongshelf components, while the across-shelf components show lower 
or statistieally insignifioant correlations, and additionally, appear more out of phase. 
These results are eonsistent with the time series plots of the mooring data (Figure 14). 

Significanee levels were ealculated following Amos and Koopmans (1963) using 
degrees of freedom estimates derived by Davis (1976). The effective degrees of freedom, 
Ncff^ is defined as 

N.ff = N/r„ 

where 

is the correlation time seale for the data. provides the degrees of freedom used in the 
eomputational proeess to determine if the individual eorrelations between the time series 
were statistieally significant with a confidence interval of 95%. Fewer degrees of 
freedom correspond to more self-eorrelated data and therefore higher significanee levels. 
Due to the lower significance levels of the across-shelf component, only the mid-depth 
eorrelation was found to be insignificant (Table 3). 

Veetor eorrelation is a teehnique that enables the average eorrelation and angular 
displaeement between a pair of two-dimensional vectors to be resolved (Kundu, 1975). 
The magnitude deseribes the overall correlation of the vector while the phase angle 
describes the average counterclockwise angle of the second vector in relation to the first 
(Kundu, 1975). Upon examination of the veetor correlations at the mooring loeations it 
beeomes evident that the veetors were well aligned and the eorrelation inereased as depth 
inereased (Table 4). 

D. SPECTRAL ANALYSIS 

Speetral analyses were eompleted to look at the time series in relation to 
frequeney. A linear detrend method was used to give the most realistie distribution of the 
spectral energy. Five overlapping pieees of the time series were used to resolve a period 
of ten days. An hourly sampling interval and a eonfidence interval of 95% were used. 
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During the principal ellipse coordinate rotation process, energy was transferred 
from the east and north components into the along and across-shelf components. This 
process resulted in an order of magnitude more energy in the alongshelf components in 
comparison to the across-shelf components (Figures 15 and 16). The alongshelf 
autospectra show that the only significant points at the 95% confidence level are the five 
and ten day periods and the tidal bands. The alongshelf autospectra for ADCP 2 also 
reveals an unambiguous inertial peak. This wind-generated peak is most visible in the 
mid-depth and near-bottom spectra. An inertial peak describes the primary frequency at 
which wind pumps energy into the ocean. The inertial period is calculated using the 
formula below for Ti where (p is latitude. 

Ti=lX (hrs) 
sin((p) 

The inertial period for this latitude is 20.16 hrs and falls between the 12 and 24 
hour tidal peaks. This inertial period is more visible at ADCP 2 than ADCP 1, indicating 
stronger wind forcing further offshore (Figures 15 and 16). The inertial oscillations at 
ADCP 1 may also be inhibited by the coastline. The diurnal tidal peak is sometimes 
broader than the semidiurnal peak at the surface with a less discernible inertial peak. Due 
to the complexities of the surface environment, some of the inertial motions in the surface 
currents is included with the broad diurnal peak. The autospectra of the near-bottom 
currents resolve the inertial peak in addition to both tidal peaks. 

The autospectra for the alongshelf and across-shelf wind were also computed 
(Figure 17). The magnitude of the across-shelf component was also smaller, by almost 
two orders of magnitude, than the alongshelf wind spectrum. The diurnal sea breeze peak 
is evident in both the across and along shore spectra. Consistent with the ADCP 
autospectra, the wind spectra also show the same significant five and ten day periods. 

E. COHERENCE AND PHASE CALCULATIONS 

The final numerical analysis applied to the wind and current time series was a test 
for coherence and phase. Similar to the autospectra, coherence and phase outputs are a 
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function of frequency. Coherence and phase were ealeulated between eaeh station (ADCP 
1 and 2) and with the respeetive aeross-shelf and alongshelf wind eomponents. The lag 
time was ealeulated using the following equation. 


f 


V 


0 


360 


X period {hrs) 


where 0 is the phase angle. 

Figure 18 depiets the eoherenee and phase between the alongshelf wind and 
alongshelf eurrent, and the aeross-shelf winds and aeross-shelf eurrent at ADCP 1 at eaeh 
depth. The straight horizontal line on the eoherenee plot indieates the 95% signifieanee 
level. While randomly signifieant, information with frequeney higher than 10'^ is 
eonsidered noise. Based on the autospeetral analysis only primary wind-band periods of 
5-10 days were eonsidered for further study. At ADCP 1 there is evidenee of surfaee 
eoherenee with a phase angle of 21° eorresponding to a lag time of 14.9 hours (Table 5). 
Signifieant eoherenees were additionally observed at both the mid-depth and near-bottom 
eurrents in the alongshelf eomponent. 

Contrary to ADCP 1, there were no signifieant wind-band periods in the aeross- 
shelf eomponents at ADCP 2 (Figure 19, Table 6). There were however signifieant 
eoherenees in all three depths of the alongshelf eomponent at ADCP 2. The mid-depth 
and near-bottom eurrents have similar phase angles and lag times, however the surfaee 
eurrent has a larger phase angle and lag time (Figures 19 and 20). 

The eoherenee and phase was also analyzed between the alongshelf wind and the 
aeross-shelf eurrent. Explained by oeean dynamies and air sea interaetions, the alongshelf 
wind forees the aeross-shelf eurrent near the surfaee. This ease was studied to diseem any 
Ekman-type flow. Eoeusing on the 5-10 day period there was signifieant eoherenee at the 
surfaee and near-bottom eurrents at ADCP 2 (Table 7). Both the 5.16 day and 10.31 day 
periods were ealeulated for the surfaee eurrents with matehing phase and lag times. The 
near-bottom eurrent was also signifieant at 10.31 days but with opposite phase. There 
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was no significant coherence for any period in the mid-depth (Figure 20). This result can 
be direetly related to Ekman dynamies. The response to the wind was offshore at the 
surfaee, onshore at the bottom, with weak flow in between. 


F. RESPONSE OF SHELF CURRENTS TO WIND EVENTS 

As seen elearly in the aireraft observations, there was submesoseale variability in 
the offshore wind distribution. During the relaxation event in the beginning of the month 
there are stronger winds north of the eape and a wind minimum at the mouth of the Bay 
(Figure 21). During the upwelling events starting on 08 August, all aireraft observations 
show an alongshelf wind maximum off Point Aho Nuevo. There is also a wind maximum 
at Point Aho Nuevo during the upwelling event from the 09-11 August (Figure 22). This 
may aeeount for the formation of an upwelling eenter there. Consistent with previous 
studies, there was a wind shadow behind the eoastal mountains in the northeast eorner of 
the Monterey Bay (Graham and Largier, 1997; Ramp et ah, 2005). 

The auto and eross-speetral analysis demonstrated some wind/eurrent interaetions 
with the alongshelf wind and the aeross-shelf eurrent. There were diurnal and semidiurnal 
peaks present as well as indieations of wind generated inertial peaks. These peaks display 
some loeal wind/sea interaetions in the time series. As usual on the eontinental shelf 
more eoherenee was found between the alongshelf wind and eurrents than between the 
aeross-shelf wind and eurrents. 

Based on the four eompleted analyses (visual, speetral, eorrelations and eoherenee 
and phase), it is evident that the eurrents appear to be foreed by both wind and non-loeal 
events sueh as eddies, meanders, and the alongshelf pressure gradient. As previously 
noted, there is some previous evidenee that the poleward pressure gradient foree 
inereases during the progression of the upwelling season (Ramp, et ah, 2005). Sinee this 
study was eompleted in August, towards the end of the upwelling season, there is reason 
to believe that the poleward pressure gradient foree is responsible for the observed 
poleward flow. This pressure gradient is driven by a eombination of two oeeanie 
phenomena. Warmer waters to the south, near Point Coneeption, result in higher eoastal 
sea levels there thus ereating a pressure gradient from the high sea surfaee levels in the 
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south to the lower sea surfaee levels in the north along the West Coast. This pressure 
gradient results in a poleward flow up the eoast. Another aspeet of this poleward flow is 
related to the stronger upwelling favorable winds located in the north. Stronger winds 
tend to set down the sea level more in the north causing a higher sea level in the south 
and a similar situation ensues, a pressure gradient is formed with resulting poleward flow. 
Both the upwelling favorable winds in the north and the warmer sea surface temperatures 
in the south contribute to the poleward flow. This is evident in the data, particularly at the 
end of the month. The wind is clearly upwelling favorable but the current flow is 
poleward. 

Given the observations it appears logical that the wind forcing would affect the 
current flow. This hypothesis could be more clearly demonstrated in an environment that 
is not interrupted by eddies and meanders and most likely influenced by the California 
Current (Figure 23). Figure 23 is an example of model outputs for the area of interest 
clearly depicting eddies and other oceanographic phenomena that could also affect the 
current flow. Previous studies have also cited that the California Current meander could 
also affect circulation patterns in the Monterey Bay region (Ramp et al., 2005; Rosenfeld 
et al., 1994). Overall, an important result of this study is that both local and remote 
forcing were active on the continental shelf near Point Ano Nuevo. 
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IV. MODEL/DATA COMPARISONS 


A. DATA PROCESSING 

Accurate modeling of the oeean eireulation is of high interest to oeeanographers 
and the Navy. Another aspeet of this projeet is to analyze eomparisons between the data 
from the ADCP moorings and three numerieal models that are designed to formulate 
eurrent predietions for the eontinental shelf and Monterey Bay. As diseussed earlier, the 
model outputs used for eomparisons are from the ROMS model, the NCOM model and 
the HOPS model. In order to build a basis for eomparison, the time series from all three 
models underwent the same low frequeney filtering and eoordinate rotation as the time 
series from the ADCPs. Onee eompleted these proeedures had similar effeets on the 
aeross and alongshelf eomponents of the model outputs. Physieally eonstrained by the 
eoast, the aeross-shelf eomponent is naturally weaker. Therefore, the models show more 
energy in the alongshelf direetion, similar to the results from the ADCP rotation. 

Oeeanographie models are able to prediet the eurrents over an area that is 
signifieantly larger than the ASAP region. The best way to organize the model data is to 
examine the temporal and spatial errors of the model relative to the observed time series. 
A uniform grid of model output was ehosen to aeeount for spatial variability surrounding 
eaeh mooring (Table 8). The eenter point of the grid eorresponds to the ADCP loeation. 
A temporal analysis was eompleted using time lags at the eenter point while a spatial 
analysis was eompleted by eomparing the ADCP data with the surrounding model grid 
points. 

In order to quantify the visual analysis it was neeessary to eompute the eorrelation 
eoeffieients between the time series of the ADCP and model data. The normalized 
eorrelation values, by definition between -1 and 1, indieate how well the model data 
matehes the ADCP data. The lag (hours) indieates the temporal mismateh between the 
two time series. For well behaved funetions, the optimal eross eorrelation value will be a 
high positive number at a small lag indieating highly eorrelated sequenees. Negative 
eorrelation values indieate that the time series are 180“ out of phase at the given lag. A 
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positive lag time indieates that the data leads the model. As the time series are temporally 
aligned the differenee between the eurves should beeome smaller both visually and 
statistically. 

The correlation coefficient calculated from the normalized cross covariance 
function is used to compare the model and the data time series. (Bendat and Piersol, 
1986). It states that dividing the cross covariance {Cxy) by the variance (cr^cr^,) of both 

sequences will yield the linear dependence between the two series x (model) and y 
(ADCP data) (Bendat and Piersol, 1986). 


pM = 






where (r) is defined as the correlation coefficient (normalized cross-covariance 

function) between x and y, C^y (x) is the cross covariance between x and y, are 

the variances of series x and y, and x is the lag. 

The root mean square (rms) error describes the mean offset (cm/s) between the 
model time series and the data time series. The root mean square error is the square root 
of the sum of the individual squared errors and is defined as 


rms error = yj[0 - 

where O is the estimated value (model time series) and O is the true value (ADCP time 
series). 

The reported rms value is the rms value associated with the maximum correlation 
at each grid point (Tables 9-16). Mathematically the lagged rms will generally be smaller 
than the original rms, assuming the lag created an improvement. 

For each grid point the maximum correlation, its lag time, and the rms error 
between the model and data were recorded for both the across-shelf and alongshelf 
components. The statistical significance of the correlations was also calculated based on 
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the degrees of freedom. The ealeulation for the degrees of freedom was derived by Davis 
(1976) and was deseribed in more detail in the wind and eurrent eorrelation statisties 
section. These parameters indicate which grid points are legitimate and should be 
considered, as well as those correlations that are insignificant. 

B. NUMERICAL MODEL ANALYSIS 

In line with the goals of the ASAP experiment, results from the three associated 
numerical models were compared to the data from the ADCP moorings. Three 
comparison techniques were used: the principal axis ellipses give estimates of the low 
frequency alignment between the models and the moored time series. Stick vector plots 
provide an intuitive approach and show differences in the magnitude and direction of the 
flow. Finally the along and across-shelf components were used in conjunction with the 
statistical calculations to determine model accuracy based on correlation coefficients, rms 
error and lag time. 

I. Principal Axis Ellipses 

Before the data were filtered all of the time series were rotated based on the 
principal axis technique. Differences in length are considered discrepancies in estimates 
of magnitude while differences in width describe inconsistencies in the energy 
distribution between the alongshelf and across-shelf components. The principal axis 
ellipses show that there was more horizontal shear in the water column than the models 
predict. The lengths of the model ellipses are generally the same at ADCP 1 and ADCP 2 
while the mooring data ellipses increased in length at ADCP 2 (Figures 24 and 25). 

The green ellipses, representing the NCOM predictions show the model 
overestimated by approximately 2 cm/s relative to the remaining models and the moored 
data at ADCP 1 (Figure 24). At ADCP 2 the NCOM model underestimated and showed a 
misalignment with respect to the mooring data (Figure 25). The NCOM ellipse (-57.2°) 
was aligned outside of the ADCP envelope (-49.5°); where as the other models were more 
evenly aligned within the ADCP envelope (Table 17). The error from the misalignment 
in the NCOM ellipse is seen in subsequent analyses. The HOPS nested (HOPSn) model 
appears to have the best alignment at both mooring locations. The ROMS model shows a 
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clear misalignment (-69.5°) at ADCP 1 and slightly improved alignment (-57.5°) at 
ADCP 2. The HOPS solo (HOPSs) model predieted signifieantly more energy in the 
aeross-shelf direetion than the other models and mooring data. 

2. Stick Vector Plots 

Stick vector plots were created to provide a concise analysis of the magnitude and 
direction of the predicted current. The currents at the mooring locations may be divided 
into three general flow patterns. There was a poleward flow from 01-09 August, an 
equatorward flow from 09-18 August and the remainder of the month (18-31 August) 
was mixed but generally poleward flow. The exception was near the surface at ADCP 2. 

In the beginning of the time series at ADCP 1, NCOM shows weak variable flow 
at the surface and poleward flow at deeper depths (Figures 26, 27 and 28). NCOM 
underestimates the first poleward flow event but badly overestimates the second two 
events (15-31 August). Additionally, at ADCP 1 NCOM predicts poleward flow patterns 
from 09-12 August with observed currents that were equatorward. The ROMS model 
shows an initial poleward flow at all three depths. These initial similarities could be 
related to a longer spin-up process. The model does not however predict the equatorward 
flow from 09-11 August. Although the flow for the second two events are slightly 
overestimated, the ROMS predicted patterns appear to be more similar to the flow 
patterns of the ADCP 1 data than NCOM. 

At ADCP 1, the predictions from the HOPSn model shows the most correlated 
flow patterns with slight overestimates in current magnitude in the surface currents and 
improvements with depth. There appears to be too much poleward flow at the beginning 
of the month (01-02 August) which could be related to initialization problems or 
problems associated with the data assimilation. 

All three depths of the HOPSs model at ADCP 1 show a strong initial poleward 
flow. The model predicts only poleward flow in the surface currents throughout the 
entire time series. The strength of the poleward flow is decreased with depth resulting in 
improved agreement near the bottom. Similar to ROMS and NCOM, HOPSs does not 
predict the equatorward flow in the middle of the time series. 
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NCOM’s current predictions for ADCP 2 were often opposite from the observed 
flow at the surface (Figures 29, 30 and 31). Flow agreement improved in the mid-depth 
and near-bottom currents. Consistent with ADCP 1 results, NCOM overestimated the 
poleward flow from 15-31 August. This may be improved in a subsequent version of the 
model (I. Shulman, personal communication, 14 February 2007). In the beginning of the 
time series ROMS surface predictions were orientated in the same direction as the 
observations although the flow patterns at the mid-depth and near-bottom were all 
poleward. After 18 August, the ROMS model was anti-correlated with the surface 
currents but the model shows slight improvements with depth. The predicted mid-depth 
and near-bottom currents show flow patterns that are opposite to and larger than the 
ADCP flow and do not predict the equatorward event. 

At ADCP 2 the current magnitude and direction of the HOPSn model appear to be 
the closest match to the mooring data. At the surface, the HOPSs model shows reasonable 
predictions. At the deeper levels the model does not predict the poleward flow. The 
magnitudes are overestimated and the flow direction is misaligned. There also appears to 
be too much energy in the across-shelf component, a result consistent with the principal 
ellipse analysis. 

3. Components and Statistical Analysis 

In terms of spatial and temporal variability, the across- and alongshore component 
plots and statistical analysis provide a more quantitative method for comparing the 
models and ADCP data. The correlation coefficients, lag times, and rms errors between 
models and observations provide an overall picture of model accuracy. The calculated 
values are presented in tables for each model at both ADCP mooring locations (Tables 9- 
16). The focus of these time series will be primarily in the alongshelf component as the 
across-shelf correlations were generally statistically insignificant. To determine 
geographical offset between the model and ADCP mooring locations, model output from 
the grid points surrounding each ADCP were compared with the data in addition to the 
grid points co-located with the ADCP. 

The comparisons at the center grid points with ADCP 1 were consistent among 

three of the models (Figures 32, 33 and 34). The ROMS and both HOPS models showed 

the same pattern of decreasing correlations with depth and smallest RMS error values in 
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the near-bottom eurrents (Tables 12-16). All depths of the alongshelf eomponent of the 
NCOM model were found to be statistieally insignilieant (Table 9). The time series for 
the eenter point at ADCP 2 show high eorrelations between HOPSn, HOPSs and ROMS 
(Figures 35, 36 and 37). The NCOM model does not appear to be as in syne with the 
ADCP data as the other models. Reeall that the prineipal ellipse for NCOM at ADCP 2 
indieated over-estimates. The ROMS and both HOPS models had high eorrelation values 
on the order of 0.6 to 0.8 (Tables 12-16). Both HOPS models however showed highest 
eorrelations at the surfaee (0.88 for HOPSn and 0.77 for HOPSs). Additionally, the 
HOPSs model shows the smallest RMS error (8.15 ern/s) at the surfaee while the smallest 
RMS error for HOPSn (6.25 ern/s) was found at the near-bottom depth. 

Better agreement ean be visually observed at both mooring loeations when the 
time lags are aeeounted for (Figures 38-43). At ADCP 1, ROMS had the smallest lags 
and the both HOPS models showed eomparable lag times while NCOM had the largest 
lag times (Figures 38, 39 and 40). The lag times for the ROMS model were all positive 
indieating that the data leads the model. NCOM and both HOPS models had negative lag 
times indieating the model leads the data. 

There was slightly more temporal offset at mooring 2 (Figures 41, 42 and 43). 
Although ROMS again had the smallest lags, the surfaee eurrents predieted a positive lag 
while the mid-depth and near-bottom eurrents showed negative lags. NCOM showed 
large lag differenees deeeasing in depth as the model predieted -39 hrs at the surfaee and 
-3 hrs in the near-bottom eurrents. HOPSn reeorded larger lag times than HOPSs and 
there was zero lag at HOPSs in the surfaee eurrents. 

Spatial aceuraey of the model predietions is observed by statistieally eomparing 
neighboring grid points to the eenter point. Three of the models showed evidenee of 
improved noweasts via spatial lags (Figures 44-49). Higher eorrelations at a grid point 
other than the eenter grid point (5) indieate that the fronts and other oeeanie features are 
misplaeed with respeet to reality. ROMS and HOPSs both had the optimal eorrelations 
aligned along the eenter of the grid (grid points 2, 5, 8) indieating no east-west spatial 
misalignment. The HOPSn model shows improvements with a shift to the southern row 
of the grid (grid points 7, 8, 9). The NCOM model predietions do not show any 
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discernible patterns in ADCP 1. At ADCP 2, the ROMS model shows spatial 
improvements to the south while the NCOM model shows improvements to the west 
(grid points 1, 4, 7) (Figures 47, 48 and 49). The HOPSn model showed optimal 
eorrelations at grid point one (northwest) at all three depths. HOPSs was the only model 
that did not show any overall spatial improvements by shifting the grid pattern indieating 
it gave the most aeeurate spatial predietions. 


C. SUMMARY OF MODEL COMPARISONS 

As a main foeus of the ASAP experiment, three numerieal models were eompared 
with observations in an effort to assess their usefulness for adaptive sampling. In general 
the NCOM model produeed the least reliable results. More loealized models with shorter 
spin-up times appeared to give better predietions. Overall differenees between the models 
and mooring data are expeeted. One of the eomplex faetors in this experiment making it 
tougher for the models is that the ADCP moorings were loeated only 6.55 km apart thus 
ereating a small seale area for model resolution. 

The HOPS models were nested on a smaller loeal seale and the highest domain 
extended only to eover the Monterey Bay/ San Franeiseo regions as well as the Point Ano 
Nuevo region. Point Ano Nuevo is reeognized to be a signifieant souree of upwelling for 
the Monterey Bay and is loeated northward of the mooring loeations (Ramp et al, 2005). 
This upwelling eenter has been shown to affeet the eireulation of the northern bay area. 
Even though all the models aeeount for this outside foreing effeet the HOPS models 
appear to prediet flow patterns more aeeurately. The HOPS models were better at 
predieting the equatorward flow event (09-18 August) in the stiek plot analysis. This 
event appeared to be the hardest for the models to deseribe. This is possibly beeause it 
resulted from an eddy rather than the wind, as indieated by the barotropie nature of the 
flow. 

In addition to resolution and domains, there are many other parameters that ean 
eause differenees between model predietions. All numerieal models are affeeted by the 
differenees in data assimilation. All the models assimilated the same glider and aireraft 
data. In addition to these inputs the models assimilated other types of information into 
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their predictions. ROMS assimilated the MBARI buoy Ml and M2 temperature and 
salinity data and NCOM assimilated the MODAS information as well. It is difficult to 
determine whether these extra assimilations were helpful or not. 

At ADCP 1 the ROMS and HOPSs models both had reasonable spatial accuracy. 
At ADCP 2, HOPSn showed the most consistent model output as it predicted the best 
improvements would occur with a northwest shift (grid point 1). This shift is completely 
opposite the southern shift the HOPSn model predicted at ADCP 1. After reviewing the 
three analysis methods, the HOPS models, in particular HOPSn, were the most accurate 
in terms of the metrics chosen. 

D. MILITARY AND OPERATIONAL APPLICATIONS 

Accurate current prediction by ocean models will benefit amphibious operations, 
special warfare operations, and mine warfare in the littoral zone. This project also 
provides information useful to oceanographers interested in understanding the combined 
effects of wind forcing and continental shelf bathymetry on ocean circulation. Model 
predictions can be improved by exploring different methods of real time data gathering, 
data assimilation, and enhanced numerical techniques. 
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APPENDIX A. FIGURES 


North 

i k 

Equatorward 

Wind 


Offshore 

Transport 



Figure 1. Illustration of coastal Ekman dynamics. Ekman transport is 90“ to the right 

of the wind stress, therefore equatorward winds off the coast of 
Californian result in offshore transport, an important factor in the seasonal 
up welling cycle. 
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Figure 2. Visual depiction of the upweiling process. (From 

http://cordellbank.noaa.gov/ environment/seasons.html, last accessed 
February 2007) 
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Figure 3. Top: HOPS two way nested modeling seheme (From Pierre Lermusiaux, 
MIT) Middle: Hierarchy of Navy Coastal Ocean Models (NCOM) in the 
Pacific (from http://www7320.nrlssc.navy.mil, last accessed May 2007) 
Bottom: Digital representation of ROMS model nesting depicting a three 
level one-way nested scheme (From http://ourocean.jpl.nasa.gov/MB06/, 
last accessed April 2007). 
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Figure 4. Illustration of Sigma coordinate method for resolving vertical levels in 
ocean models. This method follows the bottom topography and is more 
effective in coastal regions and over the continental shelf and slope, p 
represents sea surface height. Sigma coordinate systems are used by all 
three numerical models, the Harvard Ocean Prediction System (HOPS), 
Navy Coastal and Ocean Modeling (NCOM) and Regional Ocean Model 
System (ROMS) (From http;//www.oc.nps.navy.mil/ 
nom/modeling/index.html, last accessed April 2007). 
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Figure 5. Illustration of z level eoordinates. This type of vertieal grid is more 
effective in the surface layers of the water column, p represents sea 
surface height (From 

http://www.oc.nps.navy.mil/nom/modeling/index.html, last accessed April 
2007). 
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Figure 6. This representation of hybrid vertieal levels shows the sigma eoordinates 
at the surface layer and the z coordinates in the deep water and away from 
abrupt topography. (From Dr. Julie Pullen, Marine Meteorology Division, 
NRL). 
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Figure 7. ADCP moorings and MBARI M2 wind buoy locations near the Monterey 
Bay. Bathymetric contours in meters are shown by the color bar. The solid 
white line is the coastline. 
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Figure 8. Wind vectors (top), across-shelf velocity contours (middle), and 
alongshore velocity contours (bottom) at ADCP 2. 
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Figure 9. Wind vectors (top), across-shelf velocity contours (middle), and 
alongshore velocity contours (bottom) at ADCP 1. 
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Figure 10. Spatial distribution of the surfaee wind at 30 m elevation as observed by 

the NFS Twin Otter Aircraft. The magnitude is shown by both the arrows 
and the color bar. The units are in m/s. 
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Figure 11. Stick plot of the velocity vectors for (a) the wind from MBARI M2 

compared with the stick plots for the (b) surface (10 m), (c) mid-depth (34 
m) and (d) near-bottom (46 m) currents for ADCP 1. The sticks are plotted 
every six hours. Relaxations in the winds are denoted by R. Wind 
relaxation events are positive (poleward) and upwelling wind events are 
negative (equatorward). Poleward current is positive and equatorward 
current is negative. 
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ADCP 2 Winds and Currents 



Figure 12. Stick plot of the velocity vectors for (a) the wind from MBARl M2 

compared with the stick plots for the (b) surface (10 m), (c) mid-depth (34 
m) and (d) near-bottom (46 m) currents for ADCP 2. The sticks are plotted 
every six hours. Relaxations in the winds are denoted by R. Wind 
relaxation events are positive (poleward) and upwelling wind events are 
negative (equatorward). Poleward current is positive and equatorward 
current is negative. 


40 




IjitiLudr Ijililudf: 


37.40 

37.30 

37.20 

37.10 
37.00 
36.00 
36.60 
36.70 
36.60 
36.60 

36.40 

36.30 

36.20 

36.10 
36.00 
35.00 
35.80 



I 


t 




-ry. .y , . , . . . , 

-123.10-122.01>-122.70-122.50-122.30-122.10-121.00-121.70 
Longitude 

Aircraft SST (®C) and Wind (m/a) 00/03/2006 (l.l. 50m) 




20 


37.40 --I 

10 


37.30 - 

16 


37.20 - 

17 


37.10 - 

16 


37.00 - 

15 


36.00 - 

14 


36.60 - 

13 

u 

*2 

36.70 - 

12 

1 

36,60 - 

11 

36.50 

10 


36.40 - 

0 


36.30 - 

6 


36,20 - 

7 


36.10 

6 


36.00 - 

5 


35,00 - 

4 


35.00 





I ' I ' I ' I ‘ I ' I ‘ I ■ I 

-123.10- 122.00-122.70- 122.50- 122.30-122.10-121.00-121.70 
Longitude 

Aircraft SST (®C) and Wind (m/s) 08/07/2006 (l.l. 50m) 



-123.10-122.00- 122.70-122.50- 122.30-122.10- 121.00- 121.70 
Longitude 

Aircraft SST ("C) and Wind (m/s) 06/04/2006 (l.l. 50m) 







Longitude 

Aircraft SST (®C) and Wind (m/s) 08/06/2006 (l.l. 50m) 


Figure 13. Spatial distribution of the wind at 30 m altitude and sea surface 

temperatures as observed by the NFS Twin Otter Aircraft for a) August 3; 
b) August 4; c) August 7; and d) August 8, 2006. The magnitude of the 
wind is shown by the arrows and the sea surface temperatures are shown 
by the color bar. The units are in m/s. 
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Figure 14. Time series eomparing data from the two ADCPs. The blaek line is ADCP 
1 and gray line is ADCP 2. 
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Figure 15. Autospectra of the surface (10 m), mid-depth (34 m), and near-bottom (46 
m) currents for ADCP 1. Note the energy difference between the across- 
shelf (top) and alongshelf (bottom) components. 
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Figure 16. Autospectra of the surface (12.24m), mid-depth (52.24m), and near¬ 
bottom (72.24m) currents for ADCP 2. Note the energy difference 
between the across-shelf (top) and alongshelf (bottom) components. 
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Figure 18. Coherence and phase plots between wind and current for surface (10 m), 
mid-depth (34 m), and near-bottom (46 m) currents for ADCP 1. A 
positive (negative) phase means wind leads (lags) the current. The red line 
indicates a 95% significance level for the coherence (From; Amos and 
Koopmans, 1963) 
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Figure 19. Coherence and phase plots between wind and current for surface (12 m), 
mid-depth (52 m), and near-bottom (72 m) currents for ADCP 2. A 
positive (negative) phase means wind leads (lags) the current. The red line 
indicates a 95% significance level for the coherence (From; Amos and 
Koopmans, 1963). 
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Figure 20. Coherence and phase plots for alongshelf wind and across-shelf currents 
for ADCP 1 and 2. A positive (negative) phase means wind leads (lags) 
the current. The red line indicates a 95% significance level for the 
coherence (From; Amos and Koopmans, 1963) 
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Figure 21. Aircraft data during the relaxation event on 02-05 August 2006. There are 
stronger winds north of the cape and minimum winds at the mouth of the 
bay. 
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Figure 22. Upwelling event from the 09-11 August 2006. Note the maximum of 
alongshore wind off Point Ano Nuevo. 
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Figure 23. Harvard Ocean Modeling System model output showing eddies, meanders, 
and other non-local events that affect the current flow. (From 
Lermusiaux,W.G., P.J. Haley, and W.G. Leslie, 2006; ASAP; adaptive 
sampling and prediction Monterey Bay 2006, Harvard University, from 
http;//oceans.deas.harvard.edu/ASAP/index_ASAP.html, last accessed 
February 2007). 
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Figure 24. Principle variance ellipses for the numerical models and ADCP 1 mooring 
data. 
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Figure 25. Principle variance ellipses for the numerical models and ADCP 2 mooring 
data. 
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Figure 26. Stick vector plots comparing surface currents at ADCP 1 with the four 
numerical models. 
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Figure 27. Stick vector plots comparing mid-depth currents at ADCP 1 with the four 
numerical models. 
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Figure 29. Stick vector plots comparing surface currents at ADCP 2 with the four 
numerical models. 
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Figure 30. Stick vector plots comparing mid-depth currents at ADCP 2 with the four 
numerical models. 
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Figure 32. Spatial comparisons between the numerical models and ADCP time series 
for surface currents at ADCP 1. The gray line is the model time series and 
the black line is the ADCP time series. The across-shelf components form 
the left column and the alongshelf components are on the right. 
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Figure 34. Spatial comparisons between the numerical models and ADCP time series 
for near-bottom currents at ADCP 1. The gray line is the model time series 
and the black line is the ADCP time series. The across-shelf components 
form the left column and the alongshelf components are on the right. 
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Figure 35. Spatial comparisons between the numerical models and ADCP time series 
for surface currents at ADCP 2. The gray line is the model time series and 
the black line is the ADCP time series. The across-shelf components form 
the left column and the alongshelf components are on the right. 
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Figure 36. Spatial comparisons between the numerical models and ADCP time series 
for mid-depth currents at ADCP 2. The gray line is the model time series 
and the black line is the ADCP time series. The across-shelf components 
form the left column and the alongshelf components are on the right. 
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Figure 37. Spatial comparisons between the numerical models and ADCP time series 
for near-bottom currents at ADCP 2. The gray line is the model time series 
and the black line is the ADCP time series. The across-shelf components 
form the left column and the alongshelf components are on the right. 
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Figure 38. As in Figure 32, but with the series aligned better aecording to the eross- 
eorrelation results. The model output has been shifted forward/baek 
aeeording to the lag/lead estimated from the peak in the eross-eorrelation 
funetion. 
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Figure 39. As in Figure 33, but with the series aligned better aeeording to the cross- 
eorrelation results. The model output has been shifted forward/baek 
aeeording to the lag/lead estimated from the peak in the eross-eorrelation 
funetion. 
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Figure 40. As in Figure 34, but with the series aligned better aecording to the cross¬ 
correlation results. The model output has been shifted forward/back 
according to the lag/lead estimated from the peak in the cross-correlation 
function. 
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Figure 41. As in Figure 35, but with the series aligned better aecording to the cross- 
eorrelation results. The model output has been shifted forward/baek 
aeeording to the lag/lead estimated from the peak in the eross-eorrelation 
function. 
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Figure 42. As in Figure 36, but with the series aligned better aeeording to the eross- 
eorrelation results. The model output has been shifted forward/baek 
aeeording to the lag/lead estimated from the peak in the eross-eorrelation 
funetion. 
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Figure 43. As in Figure 37, but with the series aligned better aeeording to the eross- 
eorrelation results. The model output has been shifted forward/baek 
aeeording to the lag/lead estimated from the peak in the eross-eorrelation 
funetion. 
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Figure 44. Optimal space/time lagged data produeing the highest eorrelations and 

minimal rms error between the model output and the data at ADCP 1. For 
example, the alongshelf ROMS data was best at grid point 8, lagged baek 
by 6 hours. 
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Figure 45. Optimal space/time lagged data produeing the highest eorrelations and 

minimal rms error between the model output and the data at ADCP 1. For 
example, the alongshelf ROMS data was best at grid point 2, lagged baek 
by 6 hours. 
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Figure 46. Optimal space/time lagged data produeing the highest eorrelations and 

minimal rms error between the model output and the data at ADCP 1. For 
example, the alongshelf ROMS data was best at grid point 5, lagged baek 
by 10 hours. 
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Figure 47. Optimal space/time lagged data produeing the highest eorrelations and 

minimal rms error between the model output and the data at ADCP 2. For 
example, the alongshelf ROMS data was best at grid point 7, lagged by 0 
hours. 
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Figure 48. Optimal space/time lagged data produeing the highest eorrelations and 

minimal rms error between the model output and the data at ADCP 2. For 
example, the alongshelf ROMS data was best at grid point 7, lagged ahead 
by 22 hours. 
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Figure 49. Optimal space/time lagged data produeing the highest correlations and 

minimal rms error between the model output and the data at ADCP 2. For 
example, the alongshelf ROMS data was best at grid point 9, lagged back 
by 7 hours. 
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APPENDIX B. TABLES 


Summary of Key Features for the three ocean models (HOPS, ROMS 
and NCOM/ICON) used during the MB06 real-time operations 

Features 

HOPS 

Pis: Allan Robinson 
robinson@ 

pacific.deas.harvard.edu 

Pierre Lermusiaux 
pierrel@ 

pacific.harvard.edu 

ROMS 

PI: Yi Chao 

Yi.Chao 

@jpl.nasa.gov 

NCOM/ICON 

PI: Igor Shulman 

Igor.Shulman@ 

nrlssc.navy.mil 

Physical model 

ling component 

Model Domain 

2-way nested (Mon-Bay/ San 
Fr. domain and Pt AN domain) 
Four-corner points: Mon- 
Bay/San Fr. domain 

NW: 123.97W, 37.187N 

NE: 122.53W, 37.876N 

SW: 122.62W, 35.403N 

SE: 121.19W, 36.092N 

Pt AN domain 

NW: 122.84W, 37.144N 

NE: 122.33W, 37.389N 

SW: 122.44W, 36.614N 

SE: 121.93W, 36.860N 

1.5 km and 0.5 km resolutions 
Grid-point: 

Mon-Bay/San Fr.: 100x155 
PtAN: 107 X 138 

Inner domain of a 3- 
level 1 -way nested 
system 

Four-comer points: 
NW: 123.58W; 
37.432N 

NE: 122.036W; 

37.83N 

SW: 122.547W; 
35.104N 

SE: 121.048W; 

35.57N 

About 1.5 km 
resolution in 
Curvilinear-orthogonal 
grids 

Gridpoint: 179x83 

NCOM ICON (notation 
p012) model is one-way 
coupled to 9km NCOM 
California Current Model 
(NCOM CCS), which is 
one way coupled to 1/8 
degree global NCOM 
model. 

Four-comer points over 
water: 

NW:123.2W, 36.9N 
NE:122.47W,37.44N 
SW:122.5W, 35.6N 
SE:121.43W, 35.86N 
Curvilinear-orthogonal 
grid with 1 -4 km 
resolution, 81x53 grid 
points. 

NCOM frsICON 
(notation p022) 
is one way coupled to 
NCOM ICON model. 
Curvilinear-orthogonal 
grid with 

0.5 - 1.5 km res. 

114x97 grid 

Vertical levels 

30 double sigma 

32 stretched sigma 

Currently 30 sigma. 
NCOM ICON version 
with 41 sigma-Z vertical 
levels is being tested. 

Horizontal 

mixing 

Scale-selective Shapiro filter 
for numerical diffusion 

Standard second order 
diffusion 

Smagorinsky formulation 

Vertical mixing 

Vertical: Eaplacian mixing, 
with viscosities and 

K-Profile 

Parameterization 

Mellor-Yamada 2.5 
scheme 


79 




Summary of Key Features for the three ocean models (HOPS, ROMS 
and NCOM/ICON) used during the MB06 real-time operations 


diffusivities estimated using (i) 
a simple turbulent mixing- 
layer model from the surface to 
a "Ekman” depth a function 
of the local wind-stress and (ii) 
below this wind mixing layer, 
a formulation based on 
Pacanowski and Philander 
(uses the local gradient 
Richardson number Ri) 

(KPP) 


Initial conditions 

Multiple initial states and dates 
have been utilized. Presently, 
the most successful IC is an 
initialization on July 27, based 
on: 

Synoptic (Mesoscale) state: 
hydrographic data (Glider, 
NPS-SST) between Jul 22 

2006 and Jul 31 2006 (as 
available on the COOP web¬ 
site on 08-11), with a 3 day de- 
correlation time. 

Mean (Larger-scale) state: the 
above synoptic data and the Pt 
Sur CTD’s Wecoma CTD’s 
(up to Aug 8) few GTSSP 

CTD proflles for July 2006 
Initial Surface Elevation and 
Barotropic Velocities set to 
sum of the standard 
geostrophic component and 
tidal component, and of data- 
based estimates of the 
undercurrent and of the CC at 
the open boundaries of the 
Mon-Bay-San Fr domain. 

Climatology spinup to 
July 26, 2003 followed 
by COAMPS forcing 

The NCOM ICON was 
initialized on July 11th of 
2006 using fields from 
the regional model 
NCOM CCS and spun up 
with tidal forcing, with 
COAMPS 3km surface 
fluxes and with NCOM 
CCS on open boundaries. 
Starting from July 21th 
of 2006 available glider 
data are being 
assimilated. 

Surface forcing 

Surface wind 

stress 

3-km CO AMPS wind stress; 
bi-cubically interpolated to 
HOPS grid 

3-km COAMPS wind 
at 10-m computed with 
Earge and Pond drag 
coefficient 

3-km CO AMPS wind 

stress 

Air-sea heat flux 

3-km COAMPS, bi-cubically 
interpolated to HOPS grid 

3-km COAMPS 
radiative heat fluxes, 
bulk formula for 
sensible and latent heat 
fluxes 

3-km COAMPS heat 
fluxes, short-wave flux 
penetrates water column. 
There is option for 
switch to bulk 
formulations 
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Summary of Key Features for the three ocean models (HOPS, ROMS 
and NCOM/ICON) used during the MB06 real-time operations 

Air-sea fresh¬ 
water flux 

3-km COAMPS precipitation 
Evaporation derived from 
COAMPS latent heat flux 
using: Evap = 
HeatEatent/(2.5008e6- 
2.3e3*Tair); Tair from 
COAMPS. 

All fields bi-cubically 
interpolated to HOPS 

Climatology 

3-km CO AMPS (on or 
off) 

River forcing 

No 

No 

No (with the option to 
turn on) 

Tides 

Yes 

Yes; There are three 
ROMS 

nowcast/forecast 
products: (1) 1.5-km 
ROMS without tides, 
(2) 1.5-kmROMS 
with tides, and (3) 
600-m ROMS with 
tides 

Yes 

Open boundary conditions 

Open ocean 

boundary 

conditions 

0.5km Pt AN domain: BCs 
provided by 1.5km Mon- 
Bay/San Fr. domain 

1.5km Mon-Bay/San Fr. 
domain: 

Temp and Salt: combination of 
mapped data values and 
Orlanski radiation 

Internal Velocity: combination 
of data derived values and 
Orlanski radiation 

Surface pressure: combination 
of data derived values, OTIS 
tidal fields and Orlanski 
radiation 

Barotropic Momentum forcing: 
combination of data derived 
values, OTIS tidal fields and 
Orlanski radiation 

1.5-kmROMS’ BC is 
provided by 5-km 
ROMS; 

5-km ROMS’ BC is 
provided by 15-km 
ROMS; 

15-km ROMS’ BC is 
based on climatology 
and mixing Orlanski 
and Flather. 

NCOM frsICON model 
BC are provided by 
NCOM ICON model. 
NCOM ICON’BC are 
provided by NCOM CCS 
model. 

NCOMCCS’BC are 
provided by the NCOM 
global model 

Data assimilation 

Assimilation 

method 

Optimal Interpolation or Error 
Subspace Statistical Estimation 

3-dimensional 
variational method 
(3DVAR) 

NRL Ocean Data 
Assimilation System 
(NCODA) (Multivariate 
Data Assimilation) 

Assimilation 

window 

12 hours (centered on 11 and 

23 GMT, near colder and 
warmer SST times). 

6 hours at 03, 09, 15, 

21 GMT 

12 hourly updates at OOZ 
and 12Z 
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Summary of Key Features for the three ocean models (HOPS, ROMS 
and NCOM/ICON) used during the MB06 real-time operations 

Observations being assimilated in the reanalysis 

CTDs 

Pt Sur 

Wecoma 

Revelle 

Pt. Sur 

Pt. Sur 

Gliders 

SIO 

WHOI 

SIO 

WHOI 

SIO 

WHOI ASAP 

WHOI EOCO 

NPS Aircraft 

SST 

Yes 

Yes 

Yes 

AUVs 

Yes (Dorado, RU06) 

Yes 

Yes 

Moorings: 
Temperature and 
Salinity 

No 

Yes 

No 

Moorings :ADCP 

No 

No 

No 

HF radars 
(CODAR) 

No 

No 

No 

Satellite data being assimilated in the reanalysis 

SST 

No 

Yes, along with 
AMSR-E for the 15- 
km ROMS 

No in NCOM ICON 
and NCOM frsICON. 

The NCOMCCS 
assimilates TandS 
profdes and SST from 
the MODAS system. The 
MOD AS TandS profdes 
and SST are derived 
from assimilation of 
altimeter SSH data; 
MCSST (Multi-Channel 
Sea Surface 

Temperature) and 
available in-situ Tand S 
profdes. 

SSH 

No 

No in 1.5-km ROMS; 
Yes in 5-km and 15- 
km ROMS 

No in NCOM ICON, 
or NCOM frsICON. 
Satellite SSHs are 
assimilated into the 
MODAS which is 
assimilated into the 
NCOM CCS and the 
NCOM global models. 

Model Outputs 

T, S, u, V 
(including 
surface current) 

Yes, netCDF fdes are provided 
daily. Each fde contains: 
nowcast through 4 8-hour 

Yes, netCDF fdes are 
provided every 6 hours 
for nowcast and every 

Yes (with SSH), netCDF 
fdes are provided daily 
around 9AM PDT with 
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Summary of Key Features for the three ocean models (HOPS, ROMS 
and NCOM/ICON) used during the MB06 real-time operations 

provided to 
MBARI for 
analysis and 
comparisons 

forecast fields. Files are 
provided both at twice daily 
and hourly time-resolutions. 

24 hours for the 48- 
hour forecast fields 

OOZ and 12Z nowcasts 
and33-hour forecast 
fields 

LCS calculation 
and results 

Yes. 

Yes. 

No 

Web site with 
model images 
and movies 

http://oceans.deas.harvard.edu/ 
ASAP/index ASAP.html 

http://ourocean.jpl.nas 

a.gov/MB06 

http://www7320.nrlssc.n 

avy.mil/ccsnrt 


Table 1. Comparative table for the numerieal models used during the ASAP projeet (From 
http://aosn.mbari.org/ooop/, last aooessed February 2007). 
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ADCPl 

ADCP 2 

Depth 

Mean 

Std Dev 

Mean 

Std Dev 

Surface 10 (12) m 

(-0.57,2.05) 

(2.47,6.84) 

(-.68,-6.68) 

(3.41,11.70) 

Mid-depth 24 (52) m 

(0.59,4.84) 

(1.20,6.88) 

(-.93,0.57) 

(2.22,11.80) 

Near-bottom 46 (72) m 

(-0.22,4.17) 

(1.06,6.27) 

(-.80,1.69) 

(1.80,10.35) 


Table 2. Means and standard deviations from the ADCP data. Statistics are reported as 
across-shelf, alongshelf 
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Component Correlations 


Max Across-shelf 

Sgn. Level 

Max Alongshelf 

Sgn. Level 

Surface 

0.48/-3 

0.39 

0.72/3 

0.60 

Mid Depth 

0 . 32/-197 

0.38 

0.86/4 

0.62 

Near Bottom 

0.50/51 

0.39 

0.86/1 

0.58 


Table 3. Correlation values between ADCP 1 and ADCP 2. The italieized numbers are 
statistieally insignifieant and the lags (hrs) are reported as the seeond number in 
the aeross-shelf and alongshelf eomponents. 
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Vector Correlations 

Magnitude 

Orientation 

Surface 

0.7022 

5.3028 

Mid Depth 

0.8258 

-1.6107 

Near Bottom 

0.8347 

-2.617 


Table 4. Vector correlations for the cross correlation between ADCP 1 and ADCP 2. The 
magnitude describes the correlation between the two vectors while the orientation 
describes the angle between the vectors. A negative orientation indicates that 
ADCP 2 is to the right of ADCP 1. 
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ADCPl 

Across-shelf 


Period (days) 

Period (hrs) 

Coherence 

Phase (deg) 

Lag (hrs) 

Surface (10.25m) 

10.31 

247.467 

0.778 

21.73 

14.94 

Along Shelf 

Mid Depth (34.25m) 

10.31 

247.467 

0.737 

28.18 

19.37 

Near Bottom (46.25m) 

10.31 

247.467 

0.768 

25.5 

17.53 


Table 5. Significant spectral peaks in the 5-10 day band from the cross spectra between the 
across-shelf wind and the across-shelf current and the alongshelf wind and the 
alongshelf current ADCP 1. A positive (negative) phase means wind leads (lags) 
the current. No data indicates no significant coherence in the 5-10 day band. 
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ADCP 2 

Alongshelf 


Period (days) 

Period (hrs) 

Coherence 

Phase (deg) 

Lag 

(hrs) 

Surface (12.24m) 

10.31 

247.467 

0.811 

44.33 


Mid Depth (52.24m) 

10.31 

247.467 

0.841 

25.09 



5.16 

123.733 

0.863 

28.07 

9.65 

Near Bottom (72.24m) 

10.31 

247.467 

0.896 

27.48 

18.89 


5.16 

123.733 

0.78 

35.44 

12.18 


Table 6. Significant spectral peaks from the alongshelf wind and the alongshelf current at 
ADCP 2. A positive (negative) phase means wind leads (lags) the current. There 
were no significant periods in the across-shelf component. 
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ADCP 2 

Along Shelf Wind/Across Shelf Current 


Period (days) 

Period (hrs) 

Coherence 

Phase 

Lag(hrs) 

Surface 

10.31 

247.467 

0.911 

14.25 

9.80 


5.16 

123.733 

0.874 

11.97 

4.11 

Mid Depth 

NS 

NS 

NS 

NS 

NS 

Near Bottom 

10.31 

247.467 

0.804 

-139.38 

-95.81 


Table 7. Significant spectral peaks from the cross spectra between the alongshelf wind and 
the across-shelf current for both mooring locations. A positive (negative) phase 
means wind leads (lags) the current. There were no significant periods at ADCP 
1 . 
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Table 8. 



The grid point loeations for the model/data eomparisons. The ADCP is at grid 
point 5. 
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NCOM Model 
ADCP 1 

Surface Currents 10.25m 


Across-shelf Components Alongshelf Components 


Correlation 

0.2681 

0.2533 

0.433 

Correlation 

0.2403 

0.2937 

0.3138 

RMS Error 

3.2717 

3.6062 

2.4753 

RMS Error 

12.0045 

11.2595 

7.2351 

Lag time (hrs) 

29 

-54 

-45 

Lag time (hrs) 

-48 

-47 

-49 


0.3911 

0.3664 

0.4154 


0.1807 

0.1795 

0.1921 


3.1497 

2.6677 

2.8473 


13.332 

12.7174 

11.8622 


34 

35 

37 


-47 

-52 

-46 


0.5199 

0.5641 

0.49 


0.1743 

0.1316 

0.1112 


4.0476 

4.0393 

4.9431 


12.8082 

13.2805 

13.6405 


36 

39 

45 


-49 

-56 

-56 


Mid Depth 34.25m 

Across-shelf Components Alongshelf Components 


Correlation 

0.2832 

0.24 

NaN 

Correlation 

0.2007 

0.2239 

NaN 

RMS Error 

1.9993 

2.591 

NaN 

RMS Error 

10.827 

7.6353 

NaN 

Lag time (hrs) 

0 

123 

-200 

Lag time (hrs) 

-33 

-40 

-200 


0.3821 

0.2729 

0.2053 


0.1462 

0.1611 

0.143 


1.5628 

1.7069 

2.2867 


11.9803 

11.3854 

8.1157 


0 

2 

-102 


-33 

-35 

-45 


0.3177 

0.2801 

0.3557 


0.1063 

0.0997 

0.3557 


2.1686 

2.6068 

3.6774 


11.661 

12.9044 

3.6774 


0 

-101 

-100 


-35 

-33 

-100 


Near Bottom 46.25m 

Across-shelf Components Alongshelf Components 


Correlation 

0.1808 

NaN 

NaN 

Correlation 

0.2496 

NaN 

NaN 

RMS Error 

1.6629 

NaN 

NaN 

RMS Error 

8.3452 

NaN 

NaN 

Lag time (hrs) 

-97 

-200 

-200 

Lag time (hrs) 

-32 

-200 

-200 


0.3711 

0.2241 

NaN 


0.1888 

0.1783 

NaN 


1.624 

2.376 

NaN 


1.6286 

8.2956 

NaN 


0 

-35 

-200 


-26 

-35 

-200 


0.4667 

0.4083 

0.2711 


0.1202 

0.1196 

0.0922 


1.5961 

2.3711 

3.1873 


11.3017 

11.8169 

9.5545 


0 

0 

-41 


-29 

-26 

-30 


Table 9. Computed eorrelation values, RMS errors and lag times (hrs) for the NCOM 
model at ADCP 1. The shaded grid points are insignilieant and the bolded 
numbers are the highest/optimal eorrelations in the grid. 
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NCOM Model 
ADCP2 

Surface Currents 12.24m 


Across-shelf Components Alongshelf Components 


Correlation 

0.6423 

0.5742 

0.4969 

Correlation 

0.4358 

0.3659 

0.2664 

RMS Error 

4.4183 

4.002 

3.8022 

RMS Error 

13.8046 

15.1099 

16.596 

Lag time (hrs) 

1 

7 

21 

Lag time (hrs) 

-31 

-31 

-31 


0.6509 

0.6326 

0.6007 


0.4114 

0.3395 

0.2547 


4.7147 

4.3708 

3.4955 


12.9178 

14.8002 

16.3248 


6 

18 

24 


-45 

-39 

-35 


0.5959 

0.6405 

0.6785 


0.3467 

0.3467 

0.2513 


5.6864 

5.3692 

4.4585 


13.5706 

13.5568 

15.2431 


17 

25 

28 


-44 

-45 

-42 


Mid Depth 52.24m 

Across-shelf Components Alongshelf Components 


Correlation 

0.3413 

0.3532 

0.3782 

Correlation 

0.2982 

0.2539 

0.2298 

RMS Error 

2.6794 

2.37 

2.268 

RMS Error 

14.6745 

16.2325 

12.7521 

Lag time (hrs) 

78 

73 

50 

Lag time (hrs) 

-11 

-8 

-139 


0.4013 

0.3805 

0.3632 


0.276 

0.2496 

0.2438 


2.4315 

2.2224 

2.4435 


11.2143 

12.0085 

12.2776 


78 

76 

-16 


-143 

-146 

-147 


0.2554 

0.3345 

0.437 


0.3478 

0.3106 

0.2652 


2.6909 

2.7756 

2.246 


9.8979 

10.924 

11.5842 


79 

-25 

-24 


-132 

-140 

-147 


Near Bottom 72.24m 

Across-shelf Components Alongshelf Components 


Correlation 

0.2045 

0.2787 

0.2837 

Correlation 

0.4122 

0.3771 

0.3509 

RMS Error 

2.4741 

2.0634 

2.7032 

RMS Error 

13.7518 

14.875 

11.9281 

Lag time (hrs) 

54 

57 

29 

Lag time (hrs) 

-4 

-2 

-8 


0.2005 

0.1724 

0.2912 


0.3758 

0.307 

0.2379 


2.4305 

2.2128 

1.8423 


12.5527 

14.7149 

15.1506 


47 

40 

26 


-6 

-3 

-3 


0.2295 

0.1488 

0.3349 


0.2936 

0.2525 

0.1624 


2.8181 

2.3357 

1.7052 


11.9181 

13.5378 

15.175 


46 

28 

21 


-4 

0 

0 


Table 10. Computed correlation values, RMS errors and lag times (hrs) for the NCOM 
model at ADCP 2. The shaded grid points are insignificant and the bolded 
numbers are the highest/optimal correlations in the grid. 
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ROMS Model 
ADCP 1 


Surface Currents 10.25m 


Across-shelf Components Alongshelf Components 


Correlation 

0.2221 

0.3654 

0.1062 

Correlation 

0.8448 

0.8233 

0.6065 

RMS Error 

2.8892 

2.6998 

3.9701 

RMS Error 

5.0072 

7.2178 

11.6067 

Lag time (hrs) 

137 

-63 

-58 

Lag time (hrs) 

8 

-3 

-10 


0.2473 

0.3307 

0.3067 


0.8257 

0.855 

0.7213 


3.2102 

2.7299 

3.0312 


4.9397 

6.5355 

10.5216 


142 

-56 

-57 


9 

1 

-5 


0.2786 

0.2484 

0.4467 


0.7922 

0.8582 

0.7438 


3.1122 

2.696 

2.3654 


5.0231 

5.7496 

10.1405 


144 

-50 

-50 


12 

6 

-2 


Mid Depth 34.25m 
Across-shelf Components 


Alongshelf Components 


Correlation 

0.3452 

0.3429 

NaN 

Correlation 

0.6972 

0.8681 

NaN 

RMS Error 

1.9444 

1.7801 

NaN 

RMS Error 

6.7749 

3.5366 

NaN 

Lag time (hrs) 

-50 

-52 

-200 

Lag time (hrs) 

6 

5 

-200 


0.3699 

0.3958 

NaN 


0.657 

0.819 

NaN 


1.9477 

1.6331 

NaN 


7.2904 

4.9045 

NaN 


-56 

-52 

-200 


6 

4 

-200 


0.3965 

0.4006 

0.4312 


0.6178 

0.7752 

0.8527 


1.8059 

1.5593 

1.4796 


7.8495 

6.1131 

3.972 


-58 

-59 

-45 


6 

6 

5 


Near Bottom 46.25m 
Across-shelf Components 


Alongshelf Components 


Correlation 

0.5706 

NaN 

NaN 

Correlation 

0.7028 

NaN 

NaN 

RMS Error 

1.9324 

NaN 

NaN 

RMS Error 

6.6455 

NaN 

NaN 

Lag time (hrs) 

-26 

-200 

-200 

Lag time (hrs) 

12 

-200 

-200 


0.4821 

0.4372 

NaN 


0.6045 

0.7329 

NaN 


1.8358 

1.2452 

NaN 


8.1336 

4.7331 

NaN 


-31 

-30 

-200 


11 

10 

-200 


0.3596 

0.3764 

NaN 


0.5226 

0.7082 

NaN 


1.9981 

1.4151 

NaN 


9.3613 

5.6206 

NaN 


-34 

-33 

-200 


12 

12 

-200 


Table 11. Computed correlation values, RMS errors and lag times (hrs) for the ROMS 
model at ADCP 1. The shaded grid points are insignificant and the bolded 
numbers are the highest/optimal correlations in the grid. NaN (Not a Number) 
indicates no model output provided at that depth/location. 
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ROMS Model 
ADCP2 

Surface Currents 12.24m 


Across-shelf Components Alongshelf Components 


Correlation 

0.2896 

0.278 

0.227 

Correlation 

0.7936 

0.8098 

0.7522 

RMS Error 

3.5129 

3.3138 

3.1065 

RMS Error 

13.1086 

13.3681 

15.0507 

Lag time (hrs) 

-83 

-83 

-80 

Lag time (hrs) 

11 

15 

26 


0.2265 

0.3127 

0.2599 


0.8376 

0.8354 

0.8164 


3.562 

2.9007 

2.6023 


12.7267 

13.2403 

14.2621 


-160 

-160 

-161 


4 

9 

13 


0.2569 

0.3275 

0.4502 


0.8742 

0.8678 

0.8515 


3.222 

2.8674 

3.4192 


12.0086 

12.7519 

13.7439 


-161 

-161 

-4 


0 

4 

8 


Mid Depth 52.24m 

Across-shelf Components Alongshelf Components 


Correlation 

0.2563 

0.1757 

0.1174 

Correlation 

0.6748 

0.6384 

0.7675 

RMS Error 

2.9006 

2.9736 

2.7497 

RMS Error 

14.9609 

14.6554 

10.9102 

Lag time (hrs) 

101 

-157 

-41 

Lag time (hrs) 

-18 

6 

20 


0.2129 

0.2123 

0.1454 


0.7356 

0.6581 

0.7518 


2.9105 

2.9761 

2.8852 


14.5869 

15.2473 

11.8415 


97 

-62 

-26 


-22 

-11 

19 


0.3084 

0.3107 

0.2276 


0.7707 

0.7108 

0.7308 


2.7362 

2.8385 

2.9507 


14.2309 

15.2828 

13.3195 


-63 

-65 

-64 


-22 

-18 

9 


Near Bottom 72.24m 

Across-shelf Components Alongshelf Components 


Correlation 

0.22 

0.143 

NaN 

Correlation 

0.7585 

0.7203 

NaN 

RMS Error 

2.6394 

2.4234 

NaN 

RMS Error 

13.2089 

11.6544 

NaN 

Lag time (hrs) 

98 

120 

-200 

Lag time (hrs) 

-11 

0 

-200 


0.2183 

0.1633 

0.2081 


0.7876 

0.7485 

0.7742 


2.591 

2.504 

1.8553 


13.2133 

12.2142 

9.4973 


-58 

-55 

131 


-17 

-4 

6 


0.3856 

0.2857 

0.102 


0.7902 

0.7615 

0.7907 


2.3212 

2.4927 

2.0534 


13.19 

13.3766 

9.4859 


-57 

-57 

-57 


-17 

-9 

7 


Table 12. Computed eorrelation values, RMS errors and lag times (hrs) for the ROMS 
model at ADCP 2. The shaded grid points are insignilieant and the bolded 
numbers are the highest/optimal eorrelations in the grid. NaN (Not a Number) 
indieates no model output provided at that depth/loeation. 
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HOPS Solo 
ADCP 1 


Surface Currents 10.25m 


Across-shelf Components Alongshelf Components 


Correlation 

0.2059 

0.3271 

0.3709 

Correlation 

0.7877 

0.7998 

0.781 

RMS Error 

3.4818 

4.0767 

4.2207 

RMS Error 

9.3815 

9.9843 

9.8934 

Lag time (hrs) 

15 

105 

107 

Lag time (hrs) 

-18 

-18 

-20 


0.2548 

0.2958 

0.2872 


0.7945 

0.7838 

0.7638 


3.997 

4.64 

4.333 


9.9579 

10.5393 

9.9173 


103 

103 

106 


-16 

-17 

-18 


0.3417 

0.3014 

0.3117 


0.7709 

0.744 

0.7218 


4.2578 

5.148 

5.0235 


10.3139 

10.107 

8.8037 


2 

0 

-63 


-14 

-16 

-17 


Mid Depth 34.25m 

Across-shelf Components Alongshelf Components 


Correlation 

NaN 

0.2772 

0.2346 

Correlation 

NaN 

0.6173 

0.6729 

RMS Error 

NaN 

1.334 

1.5414 

RMS Error 

NaN 

5.6452 

5.5276 

Lag time (hrs) 

-200 

-102 

9 

Lag time (hrs) 

-200 

-25 

-22 


0.242 

0.1992 

0.3094 


0.6391 

0.6651 

0.677 


1.2967 

1.7392 

1.4418 


5.6217 

6.083 

6.4097 


-105 

11 

9 


-24 

-21 

-21 


0.4218 

0.2098 

0.1768 


0.6748 

0.6834 

0.683 


1.208 

1.5895 

1.8845 


6.6127 

6.4621 

6.0684 


58 

133 

124 


-22 

-20 

-18 


Near Bottom 46.25m 

Across-shelf Components Alongshelf Components 


Correlation 

NaN 

NaN 

NaN 

Correlation 

NaN 

NaN 

NaN 

RMS Error 

NaN 

NaN 

NaN 

RMS Error 

NaN 

NaN 

NaN 

Lag time (hrs) 

-200 

-200 

-200 

Lag time (hrs) 

-200 

-200 

-200 


NaN 

0.1787 

0.2876 


NaN 

0.5934 

0.6608 


NaN 

1.7706 

1.396 


NaN 

5.1045 

5.4576 


-200 

-19 

101 


-200 

-20 

-18 


0.2763 

0.3897 

0.475 


0.6298 

0.6609 

0.6575 


1.663 

1.2937 

1.3149 


5.5276 

6.4549 

6.6923 


98 

100 

101 


-21 

-18 

-16 


Table 13. Computed correlation values, RMS errors and lag times (hrs) for the HOPS Solo 
model at ADCP 1. The shaded grid points are insignificant and the bolded 
numbers are the highest/optimal correlations in the grid. NaN (Not a Number) 
indicates no model output provided at that depth/location. 
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HOPS Solo 
ADCP2 

Surface Currents 12.24m 


Across-shelf Components Alongslielf Components 


Correlation 

0.5452 

0.5454 

0.5996 

Correlation 

0.8283 

0.7946 

0.7294 

RMS Error 

4.2095 

5.4592 

6.8971 

RMS Error 

9.9418 

11.0821 

12.0466 

Lag time (hrs) 

0 

0 

0 

Lag time (hrs) 

0 

0 

27 


0.5772 

0.6049 

0.6719 


0.8106 

0.7689 

0.6721 


5.6197 

7.9138 

6.9925 


8.2623 

8.5113 

9.5232 


0 

0 

0 


-8 

0 

0 


0.669 

0.6686 

0.6978 


0.7924 

0.7236 

0.6808 


6.5606 

5.6661 

4.4766 


8.1689 

8.4868 

9.3196 


1 

4 

1 


-13 

-18 

0 


Mid Depth 52.24m 

Across-shelf Components Alongshelf Components 


Correlation 

0.3894 

0.3577 

0.2838 

Correlation 

0.6954 

0.6972 

0.7069 

RMS Error 

2.2007 

2.6369 

5.3663 

RMS Error 

9.3829 

8.5174 

9.097 

Lag time (hrs) 

63 

70 

71 

Lag time (hrs) 

-19 

-9 

-5 


0.2899 

0.3276 

0.2525 


0.7278 

0.7633 

0.7419 


2.9644 

6.083 

13.2437 


8.4102 

10.5921 

13.7329 


54 

70 

74 


-8 

-6 

-2 


0.3927 

0.4355 

0.2962 


0.7829 

0.7356 

0.7919 


3.8569 

5.8625 

24.322 


14.1991 

10.8708 

21.4509 


57 

70 

63 


-8 

-9 

-8 


Near Bottom 72.24m 

Across-shelf Components Alongshelf Components 


Correlation 

0.2819 

0.1952 

1.631 

Correlation 

0.7324 

0.7283 

0.7296 

RMS Error 

4.1452 

6.9088 

7.2566 

RMS Error 

7.942 

7.9259 

9.8922 

Lag time (hrs) 

5 

3 

70 

Lag time (hrs) 

-15 

-7 

3 


0.2786 

0.3029 

0.2083 


0.7483 

0.7553 

0.7605 


4.9049 

4.2376 

3.9341 


7.2663 

8.6404 

12.4354 


53 

10 

7 


-7 

8 

7 


0.3654 

0.3412 

0.3046 


0.781 

0.8058 

0.7399 


3.7677 

8.111 

8.0529 


13.7516 

15.5224 

13.6092 


6 

64 

85 


-3 

-3 

8 


Table 14. Computed eorrelation values, RMS errors and lag times (hrs) for the HOPS Solo 
model at ADCP 2. The shaded grid points are insignilieant and the bolded 
numbers are the highest/optimal eorrelations in the grid. 
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HOPS Nested 
ADCP 1 

Surfaee Currents 10.25m 


Aeross-shelf Components Alongshelf Components 


Correlation 

0.3291 

0.2753 

0.2767 

Correlation 

0.7136 

0.7105 

0.7083 

RMS Error 

3.5035 

3.8184 

4.7509 

RMS Error 

6.3303 

6.3654 

6.3785 

Lag time (hrs) 

-137 

-136 

2 

Lag time (hrs) 

-12 

-12 

-12 


0.3058 

0.2953 

0.3153 


0.7227 

0.7207 

0.7185 


3.6577 

4.7431 

4.9254 


6.216 

6.1842 

6.139 


-137 

2 

1 


-12 

-12 

-13 


0.3016 

0.3228 

0.3411 


0.7307 

0.7274 

0.7263 


4.6166 

4.8597 

4.9701 


6.1023 

6.0385 

5.9144 


2 

2 

1 


-12 

-13 

-14 


Mid Depth 34.25m 

Aeross-shelf Components Alongshelf Components 


Correlation 

0.5559 

0.5527 

0.5459 

Correlation 

0.6883 

0.6807 

0.6734 

RMS Error 

1.5092 

1.5456 

1.5761 

RMS Error 

5.4306 

5.4861 

5.5257 

Lag time (hrs) 

8 

8 

8 

Lag time (hrs) 

-21 

-21 

-20 


0.563 

0.5505 

0.5349 


0.6858 

0.6789 

0.675 


1.3971 

1.4304 

1.5355 


5.4713 

5.5024 

5.516 


8 

8 

9 


-21 

-20 

-20 


0.55 

0.5292 

0.5069 


0.6878 

0.6855 

0.6884 


1.3979 

1.4943 

1.5955 


5.4604 

5.4541 

5.4166 


10 

10 

11 


-20 

-20 

-20 


Near Bottom 46.25m 

Aeross-shelf Components Alongshelf Components 


Correlation 

0.397 

0.4167 

0.4475 

Correlation 

0.586 

0.6012 

0.6101 

RMS Error 

1.8324 

1.9044 

1.8699 

RMS Error 

5.7792 

5.5582 

5.4469 

Lag time (hrs) 

106 

105 

105 

Lag time (hrs) 

-21 

-20 

-20 


0.4718 

0.5011 

0.4905 


0.633 

0.6384 

0.6413 


1.5645 

1.4759 

1.4621 


5.2651 

5.2026 

5.1838 


107 

106 

104 


-21 

-21 

-21 


0.5192 

0.5043 

0.4725 


0.6578 

0.6607 

0.6692 


1.3076 

1.2976 

1.3704 


5.0833 

5.0704 

5.0161 


106 

108 

107 


-22 

-21 

-20 


Table 15. Computed eorrelation values, RMS errors and lag times (hrs) for the HOPS 

Nested model at ADCP 1. The shaded grid points are insignifieant and the bolded 
numbers are the highest/optimal eorrelations in the grid. 


97 









HOPS Nested 
ADCP2 

Surfaee Currents 12.24m 


Aeross-shelf Components Alongshelf Components 


Correlation 

0.6103 

0.6031 

0.5928 

Correlation 

0.9014 

0.8941 

0.8855 

RMS Error 

5.1774 

5.2 

5.1964 

RMS Error 

8.6438 

8.6661 

8.7924 

Lag time (hrs) 

9 

9 

9 

Lag time (hrs) 

-4 

-4 

-5 


0.6273 

0.6237 

0.6149 


0.8904 

0.8815 

0.8722 


5.4629 

5.5645 

5.5551 


8.3172 

8.3126 

8.405 


8 

8 

9 


-5 

-4 

-4 


0.6462 

0.6301 

0.6212 


0.8799 

0.8703 

0.8628 


5.465 

5.4673 

5.4451 


8.1106 

8.0469 

8.0766 


6 

7 

9 


-5 

-5 

-4 


Mid Depth 52.24m 

Aeross-shelf Components Alongshelf Components 


Correlation 

0.6322 

0.6424 

0.6181 

Correlation 

0.8574 

0.852 

0.8489 

RMS Error 

1.6753 

1.6704 

1.7361 

RMS Error 

7.7783 

7.8143 

7.8414 

Lag time (hrs) 

17 

18 

18 

Lag time (hrs) 

-21 

-21 

-21 


0.594 

0.6071 

0.5756 


0.8496 

0.8435 

0.8437 


1.812 

1.8496 

1.9116 


7.7712 

7.8281 

7.8318 


17 

18 

18 


-19 

-19 

-19 


0.5207 

0.5217 

0.5007 


0.8506 

0.8446 

0.8428 


1.9354 

1.9208 

1.9254 


7.826 

7.8386 

7.7809 


17 

17 

16 


-18 

-18 

-18 


Near Bottom 72.24m 

Aeross-shelf Components Alongshelf Components 


Correlation 

0.3864 

0.4067 

0.407 

Correlation 

0.8626 

0.8584 

0.8525 

RMS Error 

1.8494 

1.754 

1.6988 

RMS Error 

6.4165 

6.4492 

6.607 

Lag time (hrs) 

2 

1 

1 

Lag time (hrs) 

-19 

-18 

-17 


0.4207 

0.4406 

0.4218 


0.8353 

0.8283 

0.8174 


1.7244 

1.6973 

1.8189 


6.6737 

6.8247 

7.0704 


3 

5 

10 


-16 

-16 

-16 


0.4511 

0.4825 

0.5071 


0.8144 

0.8053 

0.7969 


1.9278 

2.1469 

2.1981 


7.0156 

7.2506 

7.4344 


11 

14 

17 


-13 

-13 

-14 


Table 16. Computed correlation values, RMS errors and lag times (hrs) for the HOPS 

Nested model at ADCP 2. The shaded grid points are insignificant and the bolded 
numbers are the highest/optimal correlations in the grid. 
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Table 17. 


Angles of Principle Axis Ellipses 


Mooring 1 (deg) 

Mooring 2 (deg) 

NCOM 

-61.79 

-57.24 

ROMS 

-69.51 

-57.57 

HOPSn 

-56.61 

-48.18 

HOPSs 

-63.08 

-54.21 

ADCP 

-56.28 

-49.56 


Calculated angles for the principle axis ellipses of each model and ADCP data for 
the individual mooring locations. 
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